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Universal ion trap computation on Decoherence Free Subspaces (DFS) using only two qubit 
operations is presented. The DFS is constructed for the collective dephasing model. Encoded single 
and two-qubit logical operations are implemented via the Sorensen-Molmer interaction. Alternation 
of the effective Hamiltonians for two particular phase configurations of control fields approximates 
an anisotropic exchange interaction. This is universal over suitable encodings of one logical qubit 
into three physical qubits which are also DFS under collective decoherence. 



I. INTRODUCTION 

Quantum computation is capable of cracking some 
hard problems of both classical computation and quan- 
tum physics Q. However, it meets non-trivial difficul- 
ties on its way to physical implementation Q . One per- 
sisting source of troubles is decoherence, which smears 
out the results of the quantum computation, rendering 
them useless. One promising physical implementation 
for quantum computation that possesses relatively low 
decoherence is using trapped ions [|[ Q]. In this paper, 
we present an approach for doing quantum computation 
with trapped ions that allows one to use only pairwise in- 
teraction and stay in a Decoherence Free Subspace (DFS) 
ill: H' §1 • DFS are a subspace of states whose symmetry 
makes particular error operators act like the identity on 
the subspace. 

For ion traps, an important source of errors is long 
wavelength fluctuating magnetic fields. For neighboring 
groups of ions such long wavelength magnetic fields act as 
a collective phase error. A DFS exists for the collective 
phase error which has been studied by several groups. 
1^, 1^, 1^ . Experiments have demonstrated use of this for 
protection of two-qubit [|l0[ ^ and three-qubit sys- 
tems iQ. Experimentally, Kielpinski et al. demon- 
strated that these DFS are important in the context of 
ion traps, showing an increased phase coherence when 
one logical qubit is encoded into the DFS that is con- 
structed from two physical qubits. Recently, Kielpinski 
et al. have developed a method to perform compu- 
tation on this DFS uses entangling operations between 
groups of both two and four physical qubits. A method 
to compute on a different DFS, one that protects the ion 
trap quantum computer from errors due to spontaneous 
emission, was devolped in Ref. |l5|| . 

In this paper, we show a different way to perform uni- 
versal computation on a collective dephasing DFS, using 
only entangling operations between two physical qubits. 
We explain how this could be done in linear ion traps 
using operations that have already been established. In 
Section II we briefly describe the physical system and 
then show how one can do universal computation us- 
ing certain prepared encoded states and two Sorensen- 
Molmer-type gates ffl, . In Section III , we will explain 



subspace is a DFS. In Section IV, we conclude with some 



remarks about the feasibility of such an approach com- 
pared to other schemes for universal quantum computa- 
tion in ion traps. 



II. EFFECTIVE EXCHANGE INTERACTIONS 
IN ION TRAPS 

Our system is a string of N two level ions in a lin- 
ear ion trap |§, In a linear trap the confinement in 
the radial direction is stronger than the confinement in 
the axial direction. For small displacements of the atoms 
one can treat the motion of the system as N coupled 
one-dimcnsional harmonic oscillators. Detailed calcula- 
tions of the normal modes for a large number of atoms in 
experimentally realizable traps are now widely available 
p^ . The scheme we outline below requires that each 
atom be individually addresed by a control laser. 

Coupling between atoms is brought about by coupling 
the internal state of the ions to the motional mode of the 
string of ions. A variety of gates have been suggested to 
achieve this coupling ||, |l|, |l||. For our purposes, 
the coupling gate introduced by Sorensen and Molmer 
p6[ proves most useful. 

The Sorensen-Molmer gate, described in couples 
two ions via a two photon process that virtually popu- 
lates the excited motional modes of the ion string. We 
briefly summarize their scheme here. By applying two 
lasers of opposite detuning to two ions, i and j, one cre- 
ates the following effective Hamiltonian: 



Hi 
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where fl is the on-resonant Rabi frequency of each atom, 
rj is the Lamb-Dicke parameter, A is the detuning, and cr^ 
is the Pauli spin operator x on atom i. The Lamb-Dicke 
parameter 77 is defined as the ration of the width of the 
zero point motion of the ion divided by the wavelength 
of applied field. In terms of the ion trap parameters, 

1] — \J 2NMuj J ^' where u is the frequency of the trap, 
M is the mass of the ion, A is the wavelength of the 
driving field, and 9 is the angle between the wavevector 
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By changing the phase of the incoming laser pulses by 
7r/2, one can apply the slightly modified Hamiltonian, 



Ho 



(2) 



Note that iJiand H2 commute. Therefore, by applying 
Hi followed by H2 for equal times, one can perform the 
following unitary transformation 



U = exp 



exp(-iti?'^). (3) 



is an example of an anisotropic Heisenberg exchange 
interaction (often referred to as the XY Hamiltonian). 
Previous work in our group ||, 20, 2|, 22, 23| has shown 
that such two-body exchange Hamiltonians can be uti- 
lized to perform universal quantum computation on en- 
coded qubits. 



III. ENCODED UNIVERSALITY IN ION TRAPS 



A. Universality of Anisotropic Exchange 
Interaction 



It has recently been shown by Kempe et al. [2^, |2^ 
that the anisotropic exchange (XY) interaction is suf- 
ficient to generate universal encoded single-qubit and 
two-qubit operations. In particular, these universal 
gates may be generated over encoded subspaces spanned 
by either of the following qutrit code states consisting 
of three physical qubits: C/ = {|001), |010), |100)} or 
Cii = {|110), |101), |011)}. It has been shown that the 
anisotropic exchange interaction between two physical 
qubits within either of these encodings generates the full 
sm(3) algebra, and hence also the SU{i) group, over 
the logical qutrits. Since SU{2) C SU{Z), selecting a 
qubit out of the 3 qutrit states results in universal sin- 
gle qubit operations, as long as the full sm(3) algebra 
is preserved. In addition, encoded two-qubit operations 
have been shown to be efficiently generated from the 
anisotropic exchange interaction, p^, p3| . 



B. Encoded Space is a Decoherence Free Subspace 

The encoded space we employ here is a Decoherence 
Free Subspace (DFS) with respect to collective dephas- 
ing. An operator corresponding to this error model is de- 
fined as Sz = X]fc=i ^'zi reflecting that all physical qubits 
of the computer simultaneously gather an identical ran- 
dom shift in the relative phase between their ground and 
excited level. We illustrate the DFS on a simple example 
of a logical qubit encoded into the states |01) and |10) 
of two physical qubits. Collective dephasing simultane- 
ously shifts relative phase of each of the physical two- 
level svstems as follows: 10^ 10^. 11^ -> e^'^W). The 



logical qubit stays protected against collective dephasing 
because this phase shift becomes just an unimportant 
global phase identical for both logical qubit states. An 
important source of collective dephasing error in linear 
ion traps is long wavelength fluctuations in the ambient 
magnetic fields. Work by Kielpinski et. al ||ri|, shows 
that using the two qubit collective DFS states defined to 
protect against collective dephasing leads to an increased 
fidelity limited now only by the heating of the ion mo- 
tional state. However the results of Kempe et al. |2|] 
have shown that the anisotropic exchange interaction re- 
quires at least three qubits to serve as the logical qubit, 
with two posible qutrit codes (C/ and Cn). Therefore, 
as indicated above, one needs to add a third ion in order 
to form a single logical qubit in the minimal encoding for 
the pairwise XY interaction. 

The XY Hamiltonian preserves the DFS and the 
qutrit codes. Within the DFS, the collective dephasing 
operator Sz commutes with the operator algebra gener- 
ated by the XY interaction |l] . Intuitively, all states 
from either codes are transformed identically if all their 
physical qubits are simultaneously phase shifted as de- 
scribed in the example above. The qutrit codes C/ and 
C// preserve the total number of qubits in the excited 
level |23|. Consequently, they are invariant under collec- 
tive dephasing, Sz- 

The same properties with respect to the collective de- 
phasing model are also valid for qubit encodings derived 
from qutrits |^^ . We illustrate this here on an exam- 
ple of the logical qubit encoded into the states |001) and 
|010). The relative phase between the code-words re- 
mains unchanged if all the physical qubits of value 1 are 
simultaneously shifted by an arbitrary phase kick e*"^ 
the course of the collective dephasing process 
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In the present case, we implement the XY interaction 
in steps, alternating the a* cr^ and crj^cr^ operations. Dur- 
ing this alternation, we could fall out of the DFS because 
the interactions cr^tr^ and cr* cr^ do not commute with the 
collective dephasing operator, Sz- This departure from 
the DFS can however be minimized by rapidly switching 
between the tr* cr^ and CyC^ Hamiltonians. This method 
of avoiding decoherence is similar to proposed bang-bang 
schemes To illustrate this mechanism, we make a 

simple argument here to show why this is the case. 

Assume we would like to apply the Hamiltonian H^^ 
for a time T. One could imagine first applying cr^fi^ for a 
time T then (TyCr-i^ for a time T. Or one could apply both 
operators n times, each for a time T/n: 



U 



= (e~*'"i'^-^/"e"''"«''«^/")" (4) 



Assuming only collective decoherence, a simple popula- 
tion model, described below, shows a radically improved 
decoherence time. Assume that the states outside of the 
DFS decay with some rate 7. We expect that the total 
decoherence will be proportional to the integral over time 
of 7 times the population outside of the DFS, P{t)- An 
overestimate would be to assume no decoherence and to 
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FIG. 1: A graphical demonstration of the difference in the 
population outside the DPS, p{t), obtained by using a series 
of 2 n alternating pulses, from that obtained by applying the 
a^ax interaction followed by the (Jyay interaction once (n=l). 
One can clearly see that the total integrated population out- 
side of the DFS, / P{t)dt, wiU be reduced as one increases the 
number of alternations, n. 



simply integrate over the population outside of the DFS 
during our operation. 

Let us examine the unitary operator exp[i(T^(T^t] = 
cos(i)I + i sin(t)cr^cr^. If this operator was applied to a 
state in the DFS subspace, i.e. a qubit from either C/ or 
Cji, the identity would of course preserve this codespace 
and the population that has left the codespace, p{t), 
would then be proportional to sin^(t). The integrated 
population from to t would then be proportional to 
I — I; sin(2i). Following this operator by the application 
of exp[ial.alt] will return us completely to the subspace, 
and by symmetry leads to a total integrated population 
proportional to P{t) = i — | sin(2t). 

Returning to Eq. [4], we see that by applying our op- 
erators for a time T/n, we can achieve a total integrated 
population proportional to nP{T/n) — T ~ ^sin2r/n. 
For n large compared to T, the decoherence is simply 



proportional to Therefore, one expects a quadratic 
improvement in the fidelity with the number of steps, n. 
Fig. [1]. 

To test this simple model, we calculated the evolu- 
tion of the density matrix assuming collective decoher- 
ence and a Hamiltonian that alternates n times be- 
tween cr^^crj and cr^cr^. We then evaluated the fidelity, 
/ = tryj ^/pZpb between the evolved density matrix, 
Ph, and the expected density matrix resulting from direct 
application of the XY Hamiltonian, p^. Fig. [2] shows 
that a fit of the operation error, 1 — /, on a tt pulse as 
a function of n shows a quadratic improvement, as pre- 
dicted bv our simnle nonulation model. 
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FIG. 2: The evolution of the density matrix for a single en- 
coded qubit, consisting of three physical qubits (Cj or Cii), 
was calculated using the Lindblad equation and assuming 
only collective dephasing. The resulting system density ma- 
trix pa was then compared with the expected density matrix, 
Ph, if instead one could apply the XY Hamiltonian directly. 
The plot shows operation error, 1-/, as a function of the num- 
ber of steps, n, where / is the fidelity, f = \J ^fphpa^fpl- The 
strength of the collective decoherence, 7, was set equal to the 
strength of the driving Hamiltonian, ^ , in this example. 
The fidelity was calculated after a time -^^p: ■ A fit of the 
data for n > 16 yields the expected dependence (dashed 
line) predicted by our simple population model (see text). 



IV. CONCLUSION 

We have presented a scheme for a universal ion 
trap quantum computation on Decoherence Free Sub- 
spaces using the Sorensen-Molmer two-body Hamilto- 
nian. Rapid alternation of control pulses generates 
an effective two-body interaction which approximates 
the anisotropic exchange (XY) Hamiltonian. This en- 
sures universal quantum computation over DFS pro- 
tected qubits, as a result of the anisotropic exchange 
interaction universality properties over logical qubits en- 
coded into three physical qubits |2^, [2^. The deco- 
herence time improves quadratically with the number of 
alterations, rt, for a given total pulse duration T. 

Our scheme is scalable when one assumes the array- 
based architecture proposed by Kielpinski et. al. [Q. 
In that scheme single encoded qubits are held in individ- 
ual traps. Pairs of encoded qubits are then moved into 
the same trap in order to perform two qubit operations. 
The advantage of such a scheme is that one can avoid the 
unfavorable scaling of the two qubit interaction strength 
as one puts more ions into the same trap. Furthermore, 
by encoding the ions into DFS, qubit coherences are pro- 
tected when moving the ions from trap to trap. 

There are several differences for the physical imple- 
mentation of our scheme comDared to fH. Our scheme 
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requires one extra ion per logical qubit and each ion has 
to be individually addressed. Our operations take us out 
of the DFS, although this can be mitigated by the rapid 
alternation of pulses. The alternation of pulses results in 
only a linear increase in the total number of physical op- 
erations. An advantage of our scheme is that we use only 
two qubit entangling operations, whereas the proposal of 
Kielpinski et al. [[l4| requires four-qubit entangling op- 
erations. Two-qubit operations are inherently more ro- 
bust to noise than four qubit entangling operations [ p5[ . 
One could imagine a hybrid proposal in which the qubits 
are stored in pairs of ions. Single qubit operations are 
performed using the standard scheme ^, |ll|, Then 
when the qubits are brought together to do two qubit 
gates, two additional ions are placed in the trap allowing 



one to entangle the qubits only using two physical qubit 
operations as described here. This hybrid approach may 
offer greater flexibility. 
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